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ABSTRACT 

To date, the study of high- magnification gravitational lensing effects of galaxy clusters 
has focused upon the grossly distorted, luminous arc-like features formed in massive, 
centrally condensed clusters. We investigate the formation of a different type of im- 
age, highly magnified yet undistorted, in two widely employed cluster mass density 
profiles; an isothermal sphere with a core, and a universal dark matter halo profile 
derived from numerical simulations of Navarro et al. We examine the properties of 
images of extended sources produced by these two clusters profiles, paying particular 
attention to the undistorted images. Using simple assumptions about the source and 
lens population, we estimate the relative frequency of occurrence of highly magnified, 
undistorted images and the more commonly known giant arcs. 
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1 INTRODUCTION AND MOTIVATION 

Gravitational lensing by rich galaxy clusters produces highly 
magnified, grossly distorted images of background galaxies 
(Fort & Mellier 1994, Schindler et al. 1995, Lavery 1996). 
The mechanism of production of these "giant arcs" is well 
understood; in fact, the observed arcs are commonly used 
to determine the mass distribution in the lensing clusters 
(Lavery & Henry 1988, Fort & Mellier 1994, Kneib et al. 
1995). Despite their low surface brightness, the giant arcs 
are highly visible and easily detected due to their unique 
morphology. 

However, it is possible that galaxy clusters are capa- 
ble of producing highly magnified yet undistorted images 
of background sources. Whether such images are produced 
in reasonable numbers depends on the cluster mass pro- 
files: steep density gradients produce thin arc-like images 
(Hammer 1991, Wu & Hammer 1993), whereas flat central 
cores can easily give rise to undistorted images. Because of 
their regular morphology these latter images would not be 
as prominent as arcs, and so would not be readily recognized 
and studied. 

Undistorted magnified images are the subject of the 
present paper. Our aim is two-fold. First, we ask whether 
realistic cluster profiles can produce highly magnified undis- 
torted images, or HMUs, and if so, under what conditions. 
We consider two types of clusters: an isothermal model with 
a core, and a universal dark matter cluster model derived 
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from the numerical simulations of Navarro, Frenk, & White 
(1995, 1996). The major difference between these two pro- 
files occurs in their central regions: isothermal models have a 
flat density core, whereas those drawn from numerical simu- 
lations are characterized by gradually flattening logarithmic 
slope, with decreasing radius, and are singular at the cen- 
tre. We examine the image properties produced by these 
two profiles, paying particular attention to HMUs. Second, 
we use simple assumptions about cluster and source prop- 
erties to estimate how common such HMU images might 
be. As a way to avoid some parameter dependencies, we do 
not calculate the actual frequency of occurrence of such im- 
ages, instead we estimate the ratio of the frequency of HMUs 
to giant arcs for the two different cluster mass profiles. In 
the present paper we calculate this ratio assuming circularly 
symmetric clusters, whereas real clusters exhibit substruc- 
ture which can greatly increase the clusters' cross-section 
for generating arcs (Bartelmann et al. 1995). However our 
results can be readily converted to apply to substructured 
clusters, as outlined in the Conclusions. 

Our results indicate that the statistics of HMUs are 
quite dependent on the overall form of the lensing potential, 
suggesting a diagnostic tool of cluster cores. Such a probe 
of clusters centres is welcome as current observational tech- 
niques have yet to resolve the matter distribution in the 
inner regions of clusters. Dynamical methods are generally 
inadequate in addressing the issue because of their reliance 
on assumptions (cluster is relaxed, mass follows light), as 
well as the small numbers of galaxy radial velocities (eg. Sa- 
dat 1997). X-ray methods are more promising because the 
emitting gas is in hydrostatic equilibrium with the cluster 
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gravitational potential. X-ray derived cores tend to be large; 
based on the whole sample of the EMSS clusters, Henry et 
al. (1992) estimate that the average cluster core radius is 
125/i -1 kpc. These results are in apparent contradiction with 
strong lensing analysis of clusters. Le Fevre et al. (1994), 
based on a sample of 16 brightest high-z EMSS clusters, 
conclude that arc statistics are incompatible with X-ray core 
radii, and much smaller cores are required. However, Wax- 
man & Miralda-Escude (1995), and Miralda-Escude & Babul 
(1995) pointed out that X-ray and lensing observations can 
be reconciled, and both are compatible with a singular dark 
matter potential like the universal dark matter profile, since 
multi-phase cooling flow gas in this type of a potential tends 
to be isothermal, and naturally produces cores. In any case, 
X-ray observations do not constrain cluster dark matter 
cores. That leaves us with lensing. Lensing is favored among 
other techniques because it directly probes the distribution 
of all mass regardless of its baryonic/non-baryonic nature, 
dynamical and physical state. Strong lensing observations 
place a tight upper limit on cluster cores; however it is pos- 
sible that lensing clusters with giant arcs are a biased sample 
of clusters, because steeper-than-average density profiles are 
more likely to produce thin long arcs. Consequently, cores of 
~50/i _1 kpc, or even larger, may still prove to be a common 
feature of clusters. In fact, there is tentative observational 
evidence for the existence of cluster cores: (i) C10024+1654 
is a spectacular cluster-lens with five HST-resolved images 
of a high-z galaxy. The presence of the central demagnified 
image E requires the cluster to have a finite core (Colley et 
al. 1996, Smail et al. 1996); (ii) An extended highly luminous 
z — 2.7 galaxy was recently discovered 6" from the center 
of an intervening cluster (Yee et al. 1996). Even though its 
true 'lensing status' is still uncertain, its properties can be 
explained if the lensing cluster had a core (Williams & Lewis 
1996, but see Seitz et al. 1997). Thus the existence and sizes 
of flat cluster cores is still a matter of debate. 

If found, HMU images themselves would prove useful 
for the detailed study of galaxies at high redshift. Many 
highly distorted arcs have already been studied to deduce 
the properties of the high-z sources. For example, their star 
formation rate can be derived because spectra of magnified 
sources ca n be obtained with a reasonable telescope integra- 
tion time (Ebbels et al. 1997). High-resolution HST observa- 
tions of lensed galaxies are used to infer galaxy morphologies 
(Smail et al. 1996). However, the distorted appearance of 
arcs, coupled with the complicated nature of the lens model, 
implies that lensing inversions, which are needed in order to 
derive the kinematic and structural properties of the source 
galaxy, are difficult. Ideally, one would like to observe the 
overall simpler case of highly magnified undistorted images. 



2 CLUSTER MASS PROFILE MODELS 

We consider two circularly symmetric cluster mass profiles: 
isothermal sphere with a core (ISC) and the Navarro et al. 
(1996) universal dark matter profile (NFW). 

The 2D projected mass density profile of ISC model is 
given by 



1 +p{r/r c ) 2 
' [1 + (r/r c ) 2 ] 2 " 



(1) 



We adopt p = 0.5 to obtain an isothermal sphere at large 
radii; the ISC model. However any value between and 0.5 
results in a realistic mass model; a p = represents a Plum- 
mer model [see Schneider et al. (1992) a detailed description 
of the profile] . In Equation |l|, r c is the core radius, and £0 
is the central surface mass density. 

The 3D mass density profile of the NFW model is, 



p(r) 



r(l + r/r s ) 2 ' 



(2) 



where r s is the scale radius, and p s is the mass density at 
r ~ 0.466 r s . This model is derived from N-body simulations 
of large scale structure and is found to be applicable on 
scales of 3 x 10 11 <, M 20 o/M e <J 3 x 10 15 , where M 20 o is 
the mass within a radius r2oo- Within this radius the mean 
overdensity is a factor of 200. 

We use normalised lengths, x = r/r c and x = r/r a , 
respectively, and express the surface mass density in terms 
of the critical surface density for lensing, 



k(x) = 



D„ 



4kG D ol D L 



(3) 



where Dij are the relative angular diameter distances be- 
tween the observer (o), lens (I) and source (s). 

For ISC model, the normalized surface mass density is 



k(x) = Kq 



1 + px z 

[1 +X 2 ] 2 ' 



(4) 



For NFW model, the corresponding equation was calculated 
by Bartelmann (1996), who showed that NFW model is com- 
patible with the existence of radial arcs in clusters. We re- 
produce the relevant equations here for completeness; 



k{x) = 2k s J^ 



and f(x) = 



where re s = p s r s /E c 




(5) 



(6) 



For a circularly symmetric lens, the lens equation, which 
is the relation between the source position y, the image po- 
sition x, and the deflection angle, is given by 



J*2x'K(x')dx' m {x) 



(7) 



where m(x) is the dimensionless mass interior to x, in units 
of (c 2 /AG)(D oa D ol /D ls ). 

A lensed image is distorted both radially and tangen- 
tially with respect to the centre of the lensing potential. The 
tangential distortion is the ratio of the tangential size of the 
image to that of the source; for a sufficiently small source it 
is given by the ratio of their respective distances from lens 
centre, x/y. Similarly, the radial distortion is dx/dy, if the 
source's radial extent is dy. Since magnification is just the 
ratio of the size of the image to that of the source, it is given 
by 



/'• : 



dx x 
dy y 



(8) 



The length-to-width ratio of an image whose source is cir- 
cular is, 
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L/w = 



dy x 
dx y 



(9) 



and is the most commonly used measure of image distortion. 
Note that an isothermal density profile, k(x) oc a; -1 , implies 
that dy/dx = 1 (Equation Q), i.e. images suffer no radial 
(de)magnification. Steeper profiles always result in radially 
demagnified images, while shallower profiles usually produce 
radially magnified images. 

Using Equations we can calculate all the image 
properties needed in this paper. 



3 IMAGE PROPERTIES 
3.1 ISC model 

The top panel of Figure |f| shows the relation between the 
image and source positions, for two values of «o, of an fSC 
model. Supercritical clusters, represented here by ko = 1.1 
case (solid lines), produce three images if the source impact 
parameter is smaller than the radial caustic, y r [see Eqn 
(8.42) of Schneider et al. (1992)]. Subcritical clusters, such 
as ko — 0.9 case (dashed line) always produce one image, and 
the x-y relation is one-to-one and monotonically increasing. 

Figure ^ shows the relation between image magnifica- 
tion and distortion. Two supercritical cases are shown: solid 
lines are the three images of a ko=3.0 lens, and dot-dash 
lines are for a ko = 1.1 lens. The long-dash line is the single 
image of a critical lens, Ko=1.0, while the short-dash line is 
the single image of a ko=0.9 lens. Each lens with Ko > 1 
has three branches corresponding to three images. The pri- 
mary image, which is formed at the minimum of the lensing 
potential [see Schneider (1985) and Blandford & Narayan 
(1986)], is labeled I. This image appears on the same side 
of the lens centre as the unlensed source, and is tangentially 
extended into an arc, with L/W > 1. It is always magnified 
with respect to the source. The image formed at the saddle- 
point of the lensing potential, sometimes called the 'counter 
arc', possesses reversed parity and is labeled II. The central 
image, labeled III, is formed at the maximum of the lensing 
potential, and is radially extended, i.e. has L/W < 1. Images 
II and III are on the side of the lens opposite to the location 
of the source. Note that for very centrally condensed lenses, 
i.e. those with ko > 2, image III can be demagnified, /U < 1, 
if the source is sufficiently close to the lens centre. 

The 'hidden' parameter in this plot is the source po- 
sition, y. To illustrate its influence on the location of the 
images, we have plotted the three images of a source at 
y = 0.010 (empty circles), and the three images of a source 
at y — 0.016 (solid dots), for a no — 1.1 lens. As the source 
approaches the lens centre in projection, images I and II 
get very elongated and tend to merge along the tangential 
critical line; while image III moves close towards the lens 
centre. 

Notice that the single image branch of a ko < 1 lens 
joins to and continues as the branch of the central image, 
labeled III, of the corresponding 2 — ko lens. As an example, 
the image of a ko = 0.9 lens (short-dash line) continues as 
the central image branch of a ko = 1.1 lens. This can be 
shown as follows. For image position x very close to the lens 
centre image magnification, fj, ~ (1 — ko) -2 , has the same nu- 
merical value for a kq as well as a 2 — Ko lens. The distortion 



L/W of the image located at the centre is 1, from symmetry. 
Therefore, for small x, the branches of central images of a 
/to and a 2 — ko lens meet at the same point in the log(L/W) 
vs. \og(n) diagram. The slope of both the branches at this 
point can be shown to be dlog (L/W)/dlog (fi) = —5, in- 
dependent of Kq, and p; therefore these two branches are 
continuous. 

The most visible feature in a lensing cluster is the pri- 
mary arc, since it is always well displaced from the cluster 
centre, is highly elongated and always magnified. For this 
image, the magnification is an increasing function of distor- 
tion. This is why high magnification is associated with high 
distortion of lensed images in galaxy clusters. However high 
magnification need not always imply high distortion. 

It is apparent from Figure g that subcritical clusters, 
represented here by a ko = 0.9 case, can produce highly 
magnified undistorted images^]. For example, if the source is 
located close to the centre of a Ko = 0.9 lens, its magnifica- 
tion is ~ 100, while its distortion is negligible. In fact, the 
largest L/W ratio attained by an image of any subcritical 
lens is not greater than L/W = 3, as we now show. 

The length-to- width ratio of any image in the ISC clus- 
ter model can be derived from Equations and 



L/W = 



l-«o(l + i 2 ) p - 2 (l + a; 2 [2p-l]) 

1 - K (l +X 2 Y~ 1 



(10) 



Comparing this expression for < ko < 1 to the correspond- 
ing expression for Ko = 1, one can show that the latter is 
always greater or equal to the former, for any value of x, and 
p. Therefore the maximum value of the image distortion in 
subcritical lenses of Equation [j] can be obtained from Equa- 
tion [lo] with ko = 1 and small x. This maximum value, 
L/W — 3, is independent of and p. Thus the single image 
of any subcritical lens is rather undistorted, regardless of its 
impact parameter and magnification. 

From Figure ^| it appears that supercritical clusters 
with Ko^2 can produce undistorted images with //^lO. We 
will now show that in a realistic supercritical cluster the 
probability of detecting undistorted images is negligible. To 
address this, we need to reintroduce the projected source 
position, y. The bottom panel of Figure |l| shows the distor- 
tion and image position as a function of source position for 
images of a Ko = 1.1 lens (solid lines). It may seem that im- 
ages II and III are undistorted. Image III tends to fall very 
close to the centre of the cluster, as is seen in the upper 
panel. Since a core radius, i.e. x — 1, typically corresponds 
to ~ 50/i -1 kpc, an image at x<0.1 will probably be hid- 
den within the central cluster galaxy. That leaves us with a 
small range of impact parameters in the range \og(y) ~ -2 to 
-1.8 which could produce detectable undistorted magnified 
images of type II. For typical cluster parameters that range 
translates to roughly 0.1 arcsecond and is thus smaller than a 
typical source, i.e. a galaxy several kpc across at z s ~ 0.5 — 2. 
Additionally, as the L/W ratio changes very rapidly at these 
y's, an extended image will appear distorted. 



t We 



Section 4.3 



S V.c 



our definition of magnified and undistorted later in 
, but for now we assume that such images have mag- 



nification >10, and aspect ratio <J 3. 
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Figure 1. Image position, \x\, (top panel) and length-to-width 
ratio, (bottom panel) vs. source impact parameter, y, for two ISC 
lenses: kq = 1.1 (solid lines) and kq = 0.9 (dashed line) of the ISC 
cluster model. Roman numerals indicate image type. The dotted 
lines in the bot tom panel bracket the 'undistorted' image region, 
see Section 



1.3 



3.2 NFW model 

The NFW lens model is singular at the centre for all values 
of k 3 (Equations ^|) , and so is formally supercritical. This 
profile will always produce three images if the source impact 
parameter is sufficiently small. 

Figure ^ shows the relation between the image position, 
x and source position, y. Note that the values of k s were 
picked such that the total mass within 2.5ft" 1 Mpc of these 
clusters is the s ame as that of ISC clusters with «o=X.l, and 
0.5 (see Section 4.1). The images are labeled as in Figures [j] 
and H, and have similar properties to those described in Sec- 
tion ^ll The major difference when compared to the ISC 
model is that there is no NFW case that would correspond 
to a subcritical, i.e. single-image-only ISC cluster. In fact, 
relations between image properties, y vs. x (Figure^), and 
L/W vs. /j, (Figure Q), are qualitatively the same regardless 
of the value of k s . In particular, the distortion of the pri- 
mary images of all NFW clusters is an increasing function 
of magnification; therefore giant arcs would be a common 
type of image in these clusters. 

Clusters with k s >0.3 cannot have magnified undis- 
torted images at all (Figure ^), but those with smaller val- 
ues of k s have parts of all three image branches in the HMU 
region. Since we are interested in HMUs let us derive equa- 
tions describing image properties for small 

His Cases. We will 

take advantage of the fact that for small k s models the 
triple image region moves towards small x, see top panel of 
Figure §. Equation (2.10) of Bartelmann (1996) gives func- 
tion g(x), which is proportional to the mass interior to x\ 
mix) — 4n s g(x). For small x, this simplifies to 



-(|) 2 (2 ln[|] + l), 



Figure 2. Magnification vs. distortion of images formed by an 
isothermal cluster with a core, ISC, represented by Equation [l| 
Supercritical clusters (ko = 3.0, solid lines; kq = 1.1, dot-dash 
lines) have three image branches each. Just-critical (kq = 1-0, 
long-dash line), and subcritical (kq = 0.9, short-dash line) clus- 
ters always have one image. The images are labeled as in Figure [l| 
Horizontal dotted lines at L/W = 3 and i mark the region of 
undistorted images. The three empty (solid) circles are images of 
a source at y = 0.010 (y = 0.016). Notice that for the primary im- 
ages of supercritical lenses, log (fi) is proportional to log (L/W). 
Just subcritical lenses, on the other hand, produce images of high 
magnification and small distortion. 



(12) 



and the lens equation [7] becomes, 

x 

y » x(l + 2/s»ln[— ] + k s ). 

The condition for radial critical curve is dy/dx = 0, which 
can be solved for the radius of the corresponding radial caus- 
tic in the source plane, 

Ak s e-P/a+i/P-.]). 



Ur 



(13) 



This equation shows that no matter how small k s , yr will 
remain finite, and therefore all NFW clusters can produce 3 
images. 

The corresponding expressions for magnification and 
length-to- width of the images are derived from Equations ^, 
|E| |(| and ^, assuming small x: 

A* '■ 



(1 + 3k s + 2/t s ln[|])(l + k s + 2K s ln[|] 



L/W 



(1 + 3k s +2« s ln[|]) 



(1 + Ks + 2K s ln[f ]) 



(14) 
(15) 



(11) 



These demonstrate that a line of constant L/W in Figure |j| 
intersects the primary image branches of various k s models 
at magnifications such that /i oc k~ 2 . Therefore as one goes 
to less massive clusters the magnification of primary images 
of a given fixed distortion increases. Therefore highly mag- 
nified undistorted images of infinitesimally small sources are 
possible with low-mass NFW clusters. 
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Figure 3. Image position, (top panel) and length-to-width ra- 
tio (bottom panel) vs. source impact parameter, y, for two lenses: 
K a = 0.650 (solid lines) and k s = 0.109 (dashed line) of the NFW 
cluster model. Roman numerals indicate image type. 

To determine if these images will stay undistorted for 
sources of finite size, one needs to look at how L/W changes 
with y, i.e. the bottom panel of Figure |j| For k s = 0.109 im- 
ages II and III arise only for sources at y<0.001. For typical 
r s of 300ft _1 kpc, this corresponds to ~ 0.1 arcsecond, which 
is smaller than the expected source size. Therefore images 
of type II and III would not produce HMUs. The primary 
image shows more promise; the region of /x<;10 in Figure H 
corresponds to source position range y ~ few 0.001- 0.01, 
or ~ 1 arcsecond. 



4 MODEL ASSUMPTIONS 

In the last Section we have shown that both ISC and NFW 
models can produce HMUs under certain conditions. Next, 
we need to estimate the relative frequency of occurrence 
of HMUs with the two cluster models. Assuming that the 
cluster and the sources are at fixed redshifts, this involves 
two steps: first, for a given cluster one has to integrate over 
source impact parameters, taking into account finite source 
size, and second, integrate over the distribution of clusters 
properties. 

4.1 Galaxy clusters 

We assume that galaxy clusters form a one-parameter fam- 
ily, with the core radius r c of ISC and scale radius r s of NFW 
being fixed, but possessing a range of masses. The cluster 
mass function, dN(M)/dM , is derived from two observed 
relations of X-ray selected clusters. Based on a sample of 
EMSS clusters, Henry et al. (1992) derive cluster luminosity 
function to be, dN/dLx oc L~^* x , where ax varies with red- 
shift but is approximately 0.3 between the redshifts of 0.15 
and 0.6. A relation between cluster bolometric luminosity, 



Figure 4. Magnification vs. distortion of images formed by an 
NFW cluster model, represented by Equations |i[ All clusters, 
regardless of k s have three images each if source impact param- 
eter is sufficiently small. The images are labeled as in Figure jj| 
Horizontal dotted lines at L/W = 3 and i mark the region of 
undistorted images. The three values of k s represented in this 
Figure, 0.650, 0.217, and 0.109, were picked such that these clus- 
ters have the same mass inside 2.5/i -1 Mpc as the ISC ones with 
kq=3.0, 1.0, and 0.5, respectively. 



derived from EXOSAT data, and velocity dispersion is given 
by Edge & Stewart (1991): Lx oc a// 2,9 . Combining these 
two relations with the assumption that M oc a 2 , we arrive 
at cluster mass function, 

dN(M)/dM oc M~ 39 . (16) 

The absolute normalization is irrelevant, since we will only 
be dealing with ratios. We need not assume an upper or 
lower cluster mass limits: the lower mass cutoff is effectively 
achieved because low mass clusters have a negligible lensing 
cross-section, while at the upper mass end the lensing cross- 
section increases slower than the rate at which the numbers 
of clusters decrease due to the steep slope of their mass 
function. 

We fix the core radius r c at 50/i _1 kpc, which is sub- 
stantially smaller than the derived X-ray 'core' sizes, but is 
consistent with lensing observations (Fort & Mellier 1994). 
For the NFW clusters, we fix r s at 300ft -1 kpc, which is 
close to a typical value obtained in the Navarro et al. (1996) 
simulations. 

With the physical length scale of the two cluster models 
fixed, we can now derive the correspondence between the 
one-parameter ISC and NFW models, i.e. we ask what is 
the relation between ko and n s of clusters that have the 
same mass within 2.5h~ 1 Mpc. The latter is roughly equal 
to f2oo, and corresponds to 50r c and 8.33r s , respectively. 
Using Equations [J, Q and we obtain kq ~ 4.615k s . 



© 1997 RAS, MNRAS 000, 000-000 



6 L. L. R. Williams and G. F. Lewis 



4.2 Sources 

We assume that the unlensed parent population of HMUs 
is the same as that of arcs in clusters. The half-light radii 
of the sources of arc images are almost the same as the 
observed widths of the arcs, because the cluster potential is 
not expected to distort tangential arc images in the radial 
direction. Smail et al. (1996) measure half-light radii for a 
sample of 8 HST arcs (see their Figure 5) . The average half- 
light radius is about 0.5 arcseconds, which is what we will 
adopt in the present paper. We further assume that all the 
sources are circular with a uniform surface brightness. 

At any given redshift the luminosity function (LF) of 
sources is assumed to be a power law, with the slope cor- 
responding to that of the Schechter LF, a. The value of a 
is estimated to be about 1.1 locally, but may have been 
steeper in the past, a ~ 1.5 (Ellis et al. 1995). We use both 
values below, to account for the possible range of a's de- 
pending on the type of galaxies and their evolution. The 
results are not very sensitive to a. We neglect the sources 
brighter than L*, the characteristic luminosity of Schechter 
LF, because of their small numbers. 

We need not make any further assumptions about the 
source luminosity function, as we explain below. The number 
of images of type i, where i can be HMU or arcs, for a given 
cluster characterized by ko or n s (or cluster mass M), is 
given by, 



rii oc 



V dy 



dL-Hi 



(17) 



where Hi is the Heviside step function which is 1 if the image 
selection criteria are satisfied (Section [4. 3| ) , and otherwise; 
Lum is the faintest observable luminosity. The outer integral 
is over the source impact parameter, y in the source plane. 
The integral can be written as 



KlMyWJ 



-1 



dy-Hi 



Llim 



y [Ky)] a 1 dy- H h 



(18) 



The function fi(y) is determined by the ISC or NFW model. 
The last step in the above equation is justified because the 
minimum magni ficat ion n required for a detectable image 
is large (Section 4.3), and a is always >1. Both L; im and 



L-k are functions of source redshift, but since we are only 
interested in the ratio of hhmu to riArcs, the dependency on 
these quantities cancels out. 

As the results are quite insensitive to z s we make no 
assumptions about the source redshift distribution. 
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Figure 5. Number of HMU and arc images convolved with the 
cluster mass function, as a function of cluster mass. The verti- 
cal axis has arbitrary normalization. Solid lines are for the ISC 
model, while dashed lines are for the NFW model. Source and 
lens redshifts were fixed at z s = 0.3 and Z; = 1.0; Schechter LF 
slope was assumed to be a = 1.5. 
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Table 1. The ratio of the number of HMU images to arcs in 
the ISC model, i.e. /isCi f° r a range of source, z s and lens, 
redshifts. Two values of the Schechter LF slope a were tried: 1.5, 
and 1.1 (in parenthesis). 



5 RESULTS 

Figure |E] shows the numbers of giant arcs and HMUs per 
cluster, weighted by the cluster mass function. The num- 
bers plotted along the vertical axis are proportional to 
rii ■ dN(M)/dM, where rii is given by Equation [l8|, and 
dN{M)/dM by Equation |l6| The solid and dashed lines 
represent ISC and NFW models, respectively, and arcs and 
HMUs are labeled. Here we assumed z\ = 0.3, z B — 1.0, and 



4.3 Image selection criteria 

Images are selected based on their size or morphology. To be 
selected as either a giant arc or a HMU, an image has to be 
lensed appreciably. For an arc to be detected, its L/W has 
to exceed 10, a commonly used criterium for giant luminous 
arcs. A HMU is defined as an image with central > 10, 
and L/W < 3. Additionally, the undistorted nature of the 
image is guaranteed by requiring that the change in L/W Table 2 - Same as Table § but for the NFW model, 

ratio across the image should not exceed 50%. 
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a = 1.5. Notice that for both mass profiles, HMU images 
are produc ed b y low mass clusters only, as was explained 
in Sections 3.1 and The relative numbers of HMU and 
arc images produced by ISC and NFW profiles can be un- 
derstood intuitively. The shapes of images depend on the 
mass density gradient at the location of the images (Sec- 
tion |2|). The core region of ISC is flat, hence it produces 
roughly equal radial and tangential magnifications, and is 
therefore ideal for generating HMUs. The projected density 
of the NFW profile, on the other hand, goes as -[ln| + l] at 
small x, and so does not have a flat core. It is also steeper 
than isothermal at large radii. Overall, the NFW profile is 
steeper than ISC, and thus is better than ISC at making 
giant arcs, and worse than ISC at making HMUs. 

We are interested in the ratio of HMUs to giant arcs for 
a given cluster model, i.e., 



fn 



JnuMv(M) ■ [dN(M)/dM] dM 
J nArcs(M) • [dN(M)/dM] dM ' 



(19) 



where the model is either ISC or NFW. This is just the ratio 
of the areas under the HMU and arc curves in Figure ^. For 
the parameters of the Figure, /isc = 1.16, and /nfw = 0.23. 

Table [j] shows how the /isc ratio changes with source 
and lens redshift. The numbers (in parentheses) are for the 
Schechter LF slope a of 1.5 (1.1). Table | contains the cor- 
responding values for the NFW model. It is seen from these 
two tables that the predictions are virtually independent of 
source redshift, but depend on lens redshift. The source and 
lens redshift come in only in E cr it (Equation ^), and both 
of these trends arise because of the way E c , it, which deter- 
mines cluster's lensing strength, varies with zi and z s . When 
zi is held constant, E crit attains its asymptotic value at red- 
shifts just beyond z;, and thus source redshift has very little 
effect on /. On the other hand, when z s is constant, E cr i t 
changes quite rapidly at low-moderate lens redshifts. The 
sense of the trend is also understood in terms of E cr it: when 
the lensing strength of a cluster of a certain mass increases 
due to the decrease in E cr i t , more arcs are produced, and 
hence / goes down. The weak dependency on z s is a useful 
feature because source redshift distribution is arguably the 
most uncertain of the relevant model parameters. 

The / ratios depend somewhat on the source LF slope, 
a. The dependency is not very strong because most of the 
images are magnified just by the minimum required fi (be- 
cause the cluster cross-section declines rapidly with fi), and 
so both HMUs and arcs are drawn from a rather narrow por- 
tion of the LF. The decline in cluster cross-section with \x 
is more severe for HMUs because of the additional selection 
restrictions placed on undistorted images. Therefore, when 
q is steeper and there are more faint sources, the number 
of HMU images does not increase as much as the number of 
arcs; hence / is smaller for steeper LFs. 

Dependence on cosmology is weak. Tables |l| and ^ as- 
sume flat Universe model, dominated by a cosmological con- 
stant; Q = 0.2, A = 0.8. Adopting a matter-dominated flat 
Universe, SI = 1.0, produces at most a 40% increase in /. 
Since / is the ratio of the numbers of images, it does not 
depend on the Hubble constant. 

Note that the predictions for relative numbers of HMUs 
to giant arcs are a function of chosen image selection criteria. 
For example, if minimum HMU magnification is increased 
from 10 to 30, then /isc drops by about a factor or 10 com- 



pared to the tabulated values, but virtually no HMUs are 
expected with the universal dark matter profile. 



6 CONCLUSIONS 

In this paper we considered two cluster mass density profiles, 
isothermal sphere with a core and a universal dark matter 
profile, and studied the properties of gravitationaly lensed 
images of extented sources produced by these models. We 
were particularly interested in highly magnified undistorted 
images, and have shown that both profiles can in principle 
produce such images, though when realistic cluster proper- 
ties and mass functions were folded in, ISC models proved to 
be much more efficient than NFW at making HMU images. 
Using simplified assumptions, we calculated /, the ratio of 
HMUs to giant arcs for the two models, as a function of 
source and lens parameters. To account for cluster asym- 
metries and substructure, which were not considered in this 
paper, the derived / ratios should be multiplied by the ra- 
tio of the frequency of arcs in substructured clusters to that 
in smooth symmetric clusters. The latter ratio can be ob- 
tained from numerical simulations of the type described in 
Bartelmann et al. (1995). We have shown that / is not very 
sensitive to cosmology, source luminosity function and red- 
shift distribution, and lens redshift distribution. 

In fact, the strongest dependency by far is on the clus- 
ter model: with the ISC model HMUs are on average 1.2 
times as abundant as giant arcs, whereas with the NFW 
model HMUs are only 0.2 times as common as arcs (see Ta- 
bles ^ and ^). Thus relative frequency of highly magnified 
undistorted images to giant arcs can be used to discriminate 
between isothermal clusters with cores and universal dark 
matter halo profiles, if a complete sample of clusters is ex- 
amined for HMU/arc images. Two points need to be stressed 
here. First, since HMUs are produced by clusters at the lower 
end of the mass function, M(< 2.5/i~ 1 Mpc) < 10 15 /i _1 M Q , 
a cluster sample should extend down to such masses, or 
corresponding X-ray luminosities; alternatively, model pre- 
dictions should make an allowance for a high-mass cluster 
cutoff. Second, HMUs would not be as easily detected as gi- 
ant arcs, since they would have regular morphology. If they 
are at high redshift their surface brightness will be corre- 
spondingly low. However, they should be appreciably larger 
than typical high-z galaxies, with half-light radii of about 
1.5-2 arcseconds. To find HMUs one would need to do spec- 
troscopy on extended, low surface brightness galaxies with 
regular morphology, located within 10-20 arcseconds of clus- 
ter centres, in clusters with no giant arcs. Extensive searches 
with such observational constraints have not been under- 
taken; therefore it is not surprising that highly magnified 
undistorted images, whether their number density on the 
sky is comparable to or much smaller than that of giant 
arcs, have so far escaped detection. 
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